3.24.1 \(\int (1+x)^{3/2} (a+b x) (1-x+x^2)^{3/2} \, dx\) [2301]

3.24.1.1 Optimal result
3.24.1.2 Mathematica [C] (verified)
3.24.1.3 Rubi [A] (verified)
3.24.1.4 Maple [A] (verified)
3.24.1.5 Fricas [C] (verification not implemented)
3.24.1.6 Sympy [F]
3.24.1.7 Maxima [F]
3.24.1.8 Giac [F]
3.24.1.9 Mupad [F(-1)]

3.24.1.1 Optimal result

Integrand size = 25, antiderivative size = 365 \[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\frac {54 b \sqrt {1+x} \sqrt {1-x+x^2}}{91 \left (1+\sqrt {3}+x\right )}+\frac {18 \sqrt {1+x} \sqrt {1-x+x^2} \left (91 a x+55 b x^2\right )}{5005}+\frac {2}{143} \sqrt {1+x} \sqrt {1-x+x^2} \left (13 a x+11 b x^2\right ) \left (1+x^3\right )-\frac {27 \sqrt [4]{3} \sqrt {2-\sqrt {3}} b (1+x)^{3/2} \sqrt {1-x+x^2} \sqrt {\frac {1-x+x^2}{\left (1+\sqrt {3}+x\right )^2}} E\left (\arcsin \left (\frac {1-\sqrt {3}+x}{1+\sqrt {3}+x}\right )|-7-4 \sqrt {3}\right )}{91 \sqrt {\frac {1+x}{\left (1+\sqrt {3}+x\right )^2}} \left (1+x^3\right )}+\frac {18\ 3^{3/4} \sqrt {2+\sqrt {3}} \left (91 a-55 \left (1-\sqrt {3}\right ) b\right ) (1+x)^{3/2} \sqrt {1-x+x^2} \sqrt {\frac {1-x+x^2}{\left (1+\sqrt {3}+x\right )^2}} \operatorname {EllipticF}\left (\arcsin \left (\frac {1-\sqrt {3}+x}{1+\sqrt {3}+x}\right ),-7-4 \sqrt {3}\right )}{5005 \sqrt {\frac {1+x}{\left (1+\sqrt {3}+x\right )^2}} \left (1+x^3\right )} \]

output
18/5005*(55*b*x^2+91*a*x)*(1+x)^(1/2)*(x^2-x+1)^(1/2)+2/143*(11*b*x^2+13*a 
*x)*(x^3+1)*(1+x)^(1/2)*(x^2-x+1)^(1/2)+54/91*b*(1+x)^(1/2)*(x^2-x+1)^(1/2 
)/(1+x+3^(1/2))-27/91*3^(1/4)*b*(1+x)^(3/2)*EllipticE((1+x-3^(1/2))/(1+x+3 
^(1/2)),I*3^(1/2)+2*I)*(x^2-x+1)^(1/2)*(1/2*6^(1/2)-1/2*2^(1/2))*((x^2-x+1 
)/(1+x+3^(1/2))^2)^(1/2)/(x^3+1)/((1+x)/(1+x+3^(1/2))^2)^(1/2)+18/5005*3^( 
3/4)*(1+x)^(3/2)*EllipticF((1+x-3^(1/2))/(1+x+3^(1/2)),I*3^(1/2)+2*I)*(91* 
a-55*b*(1-3^(1/2)))*(x^2-x+1)^(1/2)*(1/2*6^(1/2)+1/2*2^(1/2))*((x^2-x+1)/( 
1+x+3^(1/2))^2)^(1/2)/(x^3+1)/((1+x)/(1+x+3^(1/2))^2)^(1/2)
 
3.24.1.2 Mathematica [C] (verified)

Result contains complex when optimal does not.

Time = 32.42 (sec) , antiderivative size = 437, normalized size of antiderivative = 1.20 \[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\frac {2 x \sqrt {1+x} \sqrt {1-x+x^2} \left (91 a \left (14+5 x^3\right )+55 b x \left (16+7 x^3\right )\right )}{5005}-\frac {9 (1+x)^{3/2} \left (-\frac {660 \sqrt {-\frac {i}{3 i+\sqrt {3}}} b \left (1-x+x^2\right )}{(1+x)^2}+\frac {165 i \sqrt {2} \left (i+\sqrt {3}\right ) b \sqrt {\frac {3 i+\sqrt {3}-\frac {6 i}{1+x}}{3 i+\sqrt {3}}} \sqrt {\frac {-3 i+\sqrt {3}+\frac {6 i}{1+x}}{-3 i+\sqrt {3}}} E\left (i \text {arcsinh}\left (\frac {\sqrt {-\frac {6 i}{3 i+\sqrt {3}}}}{\sqrt {1+x}}\right )|\frac {3 i+\sqrt {3}}{3 i-\sqrt {3}}\right )}{\sqrt {1+x}}+\frac {\sqrt {2} \left (-182 i \sqrt {3} a+55 \left (3-i \sqrt {3}\right ) b\right ) \sqrt {\frac {3 i+\sqrt {3}-\frac {6 i}{1+x}}{3 i+\sqrt {3}}} \sqrt {\frac {-3 i+\sqrt {3}+\frac {6 i}{1+x}}{-3 i+\sqrt {3}}} \operatorname {EllipticF}\left (i \text {arcsinh}\left (\frac {\sqrt {-\frac {6 i}{3 i+\sqrt {3}}}}{\sqrt {1+x}}\right ),\frac {3 i+\sqrt {3}}{3 i-\sqrt {3}}\right )}{\sqrt {1+x}}\right )}{10010 \sqrt {-\frac {i}{3 i+\sqrt {3}}} \sqrt {1-x+x^2}} \]

input
Integrate[(1 + x)^(3/2)*(a + b*x)*(1 - x + x^2)^(3/2),x]
 
output
(2*x*Sqrt[1 + x]*Sqrt[1 - x + x^2]*(91*a*(14 + 5*x^3) + 55*b*x*(16 + 7*x^3 
)))/5005 - (9*(1 + x)^(3/2)*((-660*Sqrt[(-I)/(3*I + Sqrt[3])]*b*(1 - x + x 
^2))/(1 + x)^2 + ((165*I)*Sqrt[2]*(I + Sqrt[3])*b*Sqrt[(3*I + Sqrt[3] - (6 
*I)/(1 + x))/(3*I + Sqrt[3])]*Sqrt[(-3*I + Sqrt[3] + (6*I)/(1 + x))/(-3*I 
+ Sqrt[3])]*EllipticE[I*ArcSinh[Sqrt[(-6*I)/(3*I + Sqrt[3])]/Sqrt[1 + x]], 
 (3*I + Sqrt[3])/(3*I - Sqrt[3])])/Sqrt[1 + x] + (Sqrt[2]*((-182*I)*Sqrt[3 
]*a + 55*(3 - I*Sqrt[3])*b)*Sqrt[(3*I + Sqrt[3] - (6*I)/(1 + x))/(3*I + Sq 
rt[3])]*Sqrt[(-3*I + Sqrt[3] + (6*I)/(1 + x))/(-3*I + Sqrt[3])]*EllipticF[ 
I*ArcSinh[Sqrt[(-6*I)/(3*I + Sqrt[3])]/Sqrt[1 + x]], (3*I + Sqrt[3])/(3*I 
- Sqrt[3])])/Sqrt[1 + x]))/(10010*Sqrt[(-I)/(3*I + Sqrt[3])]*Sqrt[1 - x + 
x^2])
 
3.24.1.3 Rubi [A] (verified)

Time = 0.53 (sec) , antiderivative size = 337, normalized size of antiderivative = 0.92, number of steps used = 8, number of rules used = 8, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.320, Rules used = {1210, 2392, 27, 2392, 27, 2417, 759, 2416}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (x+1)^{3/2} \left (x^2-x+1\right )^{3/2} (a+b x) \, dx\)

\(\Big \downarrow \) 1210

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \int (a+b x) \left (x^3+1\right )^{3/2}dx}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 2392

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{2} \int \frac {2}{143} (13 a+11 b x) \sqrt {x^3+1}dx+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{143} \int (13 a+11 b x) \sqrt {x^3+1}dx+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 2392

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{143} \left (\frac {3}{2} \int \frac {2 (91 a+55 b x)}{35 \sqrt {x^3+1}}dx+\frac {2}{35} \sqrt {x^3+1} \left (91 a x+55 b x^2\right )\right )+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 27

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{143} \left (\frac {3}{35} \int \frac {91 a+55 b x}{\sqrt {x^3+1}}dx+\frac {2}{35} \sqrt {x^3+1} \left (91 a x+55 b x^2\right )\right )+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 2417

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{143} \left (\frac {3}{35} \left (\left (91 a-55 \left (1-\sqrt {3}\right ) b\right ) \int \frac {1}{\sqrt {x^3+1}}dx+55 b \int \frac {x-\sqrt {3}+1}{\sqrt {x^3+1}}dx\right )+\frac {2}{35} \sqrt {x^3+1} \left (91 a x+55 b x^2\right )\right )+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 759

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{143} \left (\frac {3}{35} \left (55 b \int \frac {x-\sqrt {3}+1}{\sqrt {x^3+1}}dx+\frac {2 \sqrt {2+\sqrt {3}} (x+1) \sqrt {\frac {x^2-x+1}{\left (x+\sqrt {3}+1\right )^2}} \left (91 a-55 \left (1-\sqrt {3}\right ) b\right ) \operatorname {EllipticF}\left (\arcsin \left (\frac {x-\sqrt {3}+1}{x+\sqrt {3}+1}\right ),-7-4 \sqrt {3}\right )}{\sqrt [4]{3} \sqrt {\frac {x+1}{\left (x+\sqrt {3}+1\right )^2}} \sqrt {x^3+1}}\right )+\frac {2}{35} \sqrt {x^3+1} \left (91 a x+55 b x^2\right )\right )+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

\(\Big \downarrow \) 2416

\(\displaystyle \frac {\sqrt {x+1} \sqrt {x^2-x+1} \left (\frac {9}{143} \left (\frac {3}{35} \left (\frac {2 \sqrt {2+\sqrt {3}} (x+1) \sqrt {\frac {x^2-x+1}{\left (x+\sqrt {3}+1\right )^2}} \left (91 a-55 \left (1-\sqrt {3}\right ) b\right ) \operatorname {EllipticF}\left (\arcsin \left (\frac {x-\sqrt {3}+1}{x+\sqrt {3}+1}\right ),-7-4 \sqrt {3}\right )}{\sqrt [4]{3} \sqrt {\frac {x+1}{\left (x+\sqrt {3}+1\right )^2}} \sqrt {x^3+1}}+55 b \left (\frac {2 \sqrt {x^3+1}}{x+\sqrt {3}+1}-\frac {\sqrt [4]{3} \sqrt {2-\sqrt {3}} (x+1) \sqrt {\frac {x^2-x+1}{\left (x+\sqrt {3}+1\right )^2}} E\left (\arcsin \left (\frac {x-\sqrt {3}+1}{x+\sqrt {3}+1}\right )|-7-4 \sqrt {3}\right )}{\sqrt {\frac {x+1}{\left (x+\sqrt {3}+1\right )^2}} \sqrt {x^3+1}}\right )\right )+\frac {2}{35} \sqrt {x^3+1} \left (91 a x+55 b x^2\right )\right )+\frac {2}{143} \left (x^3+1\right )^{3/2} \left (13 a x+11 b x^2\right )\right )}{\sqrt {x^3+1}}\)

input
Int[(1 + x)^(3/2)*(a + b*x)*(1 - x + x^2)^(3/2),x]
 
output
(Sqrt[1 + x]*Sqrt[1 - x + x^2]*((2*(13*a*x + 11*b*x^2)*(1 + x^3)^(3/2))/14 
3 + (9*((2*(91*a*x + 55*b*x^2)*Sqrt[1 + x^3])/35 + (3*(55*b*((2*Sqrt[1 + x 
^3])/(1 + Sqrt[3] + x) - (3^(1/4)*Sqrt[2 - Sqrt[3]]*(1 + x)*Sqrt[(1 - x + 
x^2)/(1 + Sqrt[3] + x)^2]*EllipticE[ArcSin[(1 - Sqrt[3] + x)/(1 + Sqrt[3] 
+ x)], -7 - 4*Sqrt[3]])/(Sqrt[(1 + x)/(1 + Sqrt[3] + x)^2]*Sqrt[1 + x^3])) 
 + (2*Sqrt[2 + Sqrt[3]]*(91*a - 55*(1 - Sqrt[3])*b)*(1 + x)*Sqrt[(1 - x + 
x^2)/(1 + Sqrt[3] + x)^2]*EllipticF[ArcSin[(1 - Sqrt[3] + x)/(1 + Sqrt[3] 
+ x)], -7 - 4*Sqrt[3]])/(3^(1/4)*Sqrt[(1 + x)/(1 + Sqrt[3] + x)^2]*Sqrt[1 
+ x^3])))/35))/143))/Sqrt[1 + x^3]
 

3.24.1.3.1 Defintions of rubi rules used

rule 27
Int[(a_)*(Fx_), x_Symbol] :> Simp[a   Int[Fx, x], x] /; FreeQ[a, x] &&  !Ma 
tchQ[Fx, (b_)*(Gx_) /; FreeQ[b, x]]
 

rule 759
Int[1/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = Numer[Rt[b/a, 3]], 
s = Denom[Rt[b/a, 3]]}, Simp[2*Sqrt[2 + Sqrt[3]]*(s + r*x)*(Sqrt[(s^2 - r*s 
*x + r^2*x^2)/((1 + Sqrt[3])*s + r*x)^2]/(3^(1/4)*r*Sqrt[a + b*x^3]*Sqrt[s* 
((s + r*x)/((1 + Sqrt[3])*s + r*x)^2)]))*EllipticF[ArcSin[((1 - Sqrt[3])*s 
+ r*x)/((1 + Sqrt[3])*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b}, x] & 
& PosQ[a]
 

rule 1210
Int[((d_.) + (e_.)*(x_))^(m_)*((f_.) + (g_.)*(x_))^(n_.)*((a_) + (b_.)*(x_) 
 + (c_.)*(x_)^2)^(p_), x_Symbol] :> Simp[(d + e*x)^FracPart[p]*((a + b*x + 
c*x^2)^FracPart[p]/(a*d + c*e*x^3)^FracPart[p])   Int[(f + g*x)^n*(a*d + c* 
e*x^3)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, m, p}, x] && EqQ[b*d + a*e, 
 0] && EqQ[c*d + b*e, 0] && EqQ[m, p]
 

rule 2392
Int[(Pq_)*((a_) + (b_.)*(x_)^(n_.))^(p_), x_Symbol] :> Module[{q = Expon[Pq 
, x], i}, Simp[(a + b*x^n)^p*Sum[Coeff[Pq, x, i]*(x^(i + 1)/(n*p + i + 1)), 
 {i, 0, q}], x] + Simp[a*n*p   Int[(a + b*x^n)^(p - 1)*Sum[Coeff[Pq, x, i]* 
(x^i/(n*p + i + 1)), {i, 0, q}], x], x]] /; FreeQ[{a, b}, x] && PolyQ[Pq, x 
] && IGtQ[(n - 1)/2, 0] && GtQ[p, 0]
 

rule 2416
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Simplify[(1 - Sqrt[3])*(d/c)]], s = Denom[Simplify[(1 - Sqrt[3])*(d/c) 
]]}, Simp[2*d*s^3*(Sqrt[a + b*x^3]/(a*r^2*((1 + Sqrt[3])*s + r*x))), x] - S 
imp[3^(1/4)*Sqrt[2 - Sqrt[3]]*d*s*(s + r*x)*(Sqrt[(s^2 - r*s*x + r^2*x^2)/( 
(1 + Sqrt[3])*s + r*x)^2]/(r^2*Sqrt[a + b*x^3]*Sqrt[s*((s + r*x)/((1 + Sqrt 
[3])*s + r*x)^2)]))*EllipticE[ArcSin[((1 - Sqrt[3])*s + r*x)/((1 + Sqrt[3]) 
*s + r*x)], -7 - 4*Sqrt[3]], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && Eq 
Q[b*c^3 - 2*(5 - 3*Sqrt[3])*a*d^3, 0]
 

rule 2417
Int[((c_) + (d_.)*(x_))/Sqrt[(a_) + (b_.)*(x_)^3], x_Symbol] :> With[{r = N 
umer[Rt[b/a, 3]], s = Denom[Rt[b/a, 3]]}, Simp[(c*r - (1 - Sqrt[3])*d*s)/r 
  Int[1/Sqrt[a + b*x^3], x], x] + Simp[d/r   Int[((1 - Sqrt[3])*s + r*x)/Sq 
rt[a + b*x^3], x], x]] /; FreeQ[{a, b, c, d}, x] && PosQ[a] && NeQ[b*c^3 - 
2*(5 - 3*Sqrt[3])*a*d^3, 0]
 
3.24.1.4 Maple [A] (verified)

Time = 0.72 (sec) , antiderivative size = 360, normalized size of antiderivative = 0.99

method result size
risch \(\frac {2 x \left (385 b \,x^{4}+455 a \,x^{3}+880 b x +1274 a \right ) \sqrt {1+x}\, \sqrt {x^{2}-x +1}}{5005}+\frac {\left (\frac {54 a \left (\frac {3}{2}-\frac {i \sqrt {3}}{2}\right ) \sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}-\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}+\frac {i \sqrt {3}}{2}}}\, F\left (\sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}, \sqrt {\frac {-\frac {3}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\right )}{55 \sqrt {x^{3}+1}}+\frac {54 b \left (\frac {3}{2}-\frac {i \sqrt {3}}{2}\right ) \sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}-\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}+\frac {i \sqrt {3}}{2}}}\, \left (\left (-\frac {3}{2}-\frac {i \sqrt {3}}{2}\right ) E\left (\sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}, \sqrt {\frac {-\frac {3}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\right )+\left (\frac {1}{2}+\frac {i \sqrt {3}}{2}\right ) F\left (\sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}, \sqrt {\frac {-\frac {3}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\right )\right )}{91 \sqrt {x^{3}+1}}\right ) \sqrt {\left (1+x \right ) \left (x^{2}-x +1\right )}}{\sqrt {1+x}\, \sqrt {x^{2}-x +1}}\) \(360\)
elliptic \(\frac {\sqrt {\left (1+x \right ) \left (x^{2}-x +1\right )}\, \sqrt {1+x}\, \sqrt {x^{2}-x +1}\, \left (\frac {2 b \,x^{5} \sqrt {x^{3}+1}}{13}+\frac {2 a \,x^{4} \sqrt {x^{3}+1}}{11}+\frac {32 b \,x^{2} \sqrt {x^{3}+1}}{91}+\frac {28 a x \sqrt {x^{3}+1}}{55}+\frac {54 a \left (\frac {3}{2}-\frac {i \sqrt {3}}{2}\right ) \sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}-\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}+\frac {i \sqrt {3}}{2}}}\, F\left (\sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}, \sqrt {\frac {-\frac {3}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\right )}{55 \sqrt {x^{3}+1}}+\frac {54 b \left (\frac {3}{2}-\frac {i \sqrt {3}}{2}\right ) \sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}-\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\, \sqrt {\frac {x -\frac {1}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}+\frac {i \sqrt {3}}{2}}}\, \left (\left (-\frac {3}{2}-\frac {i \sqrt {3}}{2}\right ) E\left (\sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}, \sqrt {\frac {-\frac {3}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\right )+\left (\frac {1}{2}+\frac {i \sqrt {3}}{2}\right ) F\left (\sqrt {\frac {1+x}{\frac {3}{2}-\frac {i \sqrt {3}}{2}}}, \sqrt {\frac {-\frac {3}{2}+\frac {i \sqrt {3}}{2}}{-\frac {3}{2}-\frac {i \sqrt {3}}{2}}}\right )\right )}{91 \sqrt {x^{3}+1}}\right )}{x^{3}+1}\) \(378\)
default \(-\frac {\sqrt {1+x}\, \sqrt {x^{2}-x +1}\, \left (-770 b \,x^{8}-910 a \,x^{7}+2457 i \sqrt {3}\, \sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}\, \sqrt {\frac {i \sqrt {3}-2 x +1}{i \sqrt {3}+3}}\, \sqrt {\frac {i \sqrt {3}+2 x -1}{-3+i \sqrt {3}}}\, F\left (\sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}, \sqrt {-\frac {-3+i \sqrt {3}}{i \sqrt {3}+3}}\right ) a -1485 i \sqrt {3}\, \sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}\, \sqrt {\frac {i \sqrt {3}-2 x +1}{i \sqrt {3}+3}}\, \sqrt {\frac {i \sqrt {3}+2 x -1}{-3+i \sqrt {3}}}\, F\left (\sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}, \sqrt {-\frac {-3+i \sqrt {3}}{i \sqrt {3}+3}}\right ) b -2530 b \,x^{5}-7371 \sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}\, \sqrt {\frac {i \sqrt {3}-2 x +1}{i \sqrt {3}+3}}\, \sqrt {\frac {i \sqrt {3}+2 x -1}{-3+i \sqrt {3}}}\, F\left (\sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}, \sqrt {-\frac {-3+i \sqrt {3}}{i \sqrt {3}+3}}\right ) a -4455 \sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}\, \sqrt {\frac {i \sqrt {3}-2 x +1}{i \sqrt {3}+3}}\, \sqrt {\frac {i \sqrt {3}+2 x -1}{-3+i \sqrt {3}}}\, F\left (\sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}, \sqrt {-\frac {-3+i \sqrt {3}}{i \sqrt {3}+3}}\right ) b +8910 \sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}\, \sqrt {\frac {i \sqrt {3}-2 x +1}{i \sqrt {3}+3}}\, \sqrt {\frac {i \sqrt {3}+2 x -1}{-3+i \sqrt {3}}}\, E\left (\sqrt {-\frac {2 \left (1+x \right )}{-3+i \sqrt {3}}}, \sqrt {-\frac {-3+i \sqrt {3}}{i \sqrt {3}+3}}\right ) b -3458 a \,x^{4}-1760 b \,x^{2}-2548 a x \right )}{5005 \left (x^{3}+1\right )}\) \(608\)

input
int((1+x)^(3/2)*(b*x+a)*(x^2-x+1)^(3/2),x,method=_RETURNVERBOSE)
 
output
2/5005*x*(385*b*x^4+455*a*x^3+880*b*x+1274*a)*(1+x)^(1/2)*(x^2-x+1)^(1/2)+ 
(54/55*a*(3/2-1/2*I*3^(1/2))*((1+x)/(3/2-1/2*I*3^(1/2)))^(1/2)*((x-1/2-1/2 
*I*3^(1/2))/(-3/2-1/2*I*3^(1/2)))^(1/2)*((x-1/2+1/2*I*3^(1/2))/(-3/2+1/2*I 
*3^(1/2)))^(1/2)/(x^3+1)^(1/2)*EllipticF(((1+x)/(3/2-1/2*I*3^(1/2)))^(1/2) 
,((-3/2+1/2*I*3^(1/2))/(-3/2-1/2*I*3^(1/2)))^(1/2))+54/91*b*(3/2-1/2*I*3^( 
1/2))*((1+x)/(3/2-1/2*I*3^(1/2)))^(1/2)*((x-1/2-1/2*I*3^(1/2))/(-3/2-1/2*I 
*3^(1/2)))^(1/2)*((x-1/2+1/2*I*3^(1/2))/(-3/2+1/2*I*3^(1/2)))^(1/2)/(x^3+1 
)^(1/2)*((-3/2-1/2*I*3^(1/2))*EllipticE(((1+x)/(3/2-1/2*I*3^(1/2)))^(1/2), 
((-3/2+1/2*I*3^(1/2))/(-3/2-1/2*I*3^(1/2)))^(1/2))+(1/2+1/2*I*3^(1/2))*Ell 
ipticF(((1+x)/(3/2-1/2*I*3^(1/2)))^(1/2),((-3/2+1/2*I*3^(1/2))/(-3/2-1/2*I 
*3^(1/2)))^(1/2))))*((1+x)*(x^2-x+1))^(1/2)/(1+x)^(1/2)/(x^2-x+1)^(1/2)
 
3.24.1.5 Fricas [C] (verification not implemented)

Result contains higher order function than in optimal. Order 9 vs. order 4.

Time = 0.08 (sec) , antiderivative size = 58, normalized size of antiderivative = 0.16 \[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\frac {2}{5005} \, {\left (385 \, b x^{5} + 455 \, a x^{4} + 880 \, b x^{2} + 1274 \, a x\right )} \sqrt {x^{2} - x + 1} \sqrt {x + 1} + \frac {54}{55} \, a {\rm weierstrassPInverse}\left (0, -4, x\right ) - \frac {54}{91} \, b {\rm weierstrassZeta}\left (0, -4, {\rm weierstrassPInverse}\left (0, -4, x\right )\right ) \]

input
integrate((1+x)^(3/2)*(b*x+a)*(x^2-x+1)^(3/2),x, algorithm="fricas")
 
output
2/5005*(385*b*x^5 + 455*a*x^4 + 880*b*x^2 + 1274*a*x)*sqrt(x^2 - x + 1)*sq 
rt(x + 1) + 54/55*a*weierstrassPInverse(0, -4, x) - 54/91*b*weierstrassZet 
a(0, -4, weierstrassPInverse(0, -4, x))
 
3.24.1.6 Sympy [F]

\[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\int \left (a + b x\right ) \left (x + 1\right )^{\frac {3}{2}} \left (x^{2} - x + 1\right )^{\frac {3}{2}}\, dx \]

input
integrate((1+x)**(3/2)*(b*x+a)*(x**2-x+1)**(3/2),x)
 
output
Integral((a + b*x)*(x + 1)**(3/2)*(x**2 - x + 1)**(3/2), x)
 
3.24.1.7 Maxima [F]

\[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\int { {\left (b x + a\right )} {\left (x^{2} - x + 1\right )}^{\frac {3}{2}} {\left (x + 1\right )}^{\frac {3}{2}} \,d x } \]

input
integrate((1+x)^(3/2)*(b*x+a)*(x^2-x+1)^(3/2),x, algorithm="maxima")
 
output
integrate((b*x + a)*(x^2 - x + 1)^(3/2)*(x + 1)^(3/2), x)
 
3.24.1.8 Giac [F]

\[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\int { {\left (b x + a\right )} {\left (x^{2} - x + 1\right )}^{\frac {3}{2}} {\left (x + 1\right )}^{\frac {3}{2}} \,d x } \]

input
integrate((1+x)^(3/2)*(b*x+a)*(x^2-x+1)^(3/2),x, algorithm="giac")
 
output
integrate((b*x + a)*(x^2 - x + 1)^(3/2)*(x + 1)^(3/2), x)
 
3.24.1.9 Mupad [F(-1)]

Timed out. \[ \int (1+x)^{3/2} (a+b x) \left (1-x+x^2\right )^{3/2} \, dx=\int {\left (x+1\right )}^{3/2}\,\left (a+b\,x\right )\,{\left (x^2-x+1\right )}^{3/2} \,d x \]

input
int((x + 1)^(3/2)*(a + b*x)*(x^2 - x + 1)^(3/2),x)
 
output
int((x + 1)^(3/2)*(a + b*x)*(x^2 - x + 1)^(3/2), x)